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bstract 

The thermomagnetic instability of the critical state in superconductors is analysed with account 
of the dissipation and dispersion. The possibility is demonstrated of the existance of a nonlinear 
shok wave describing the final stage of the instability evolution in a superconductor. The structures 
possess a finite-amplitude and propagate at a constant velocity. The apperance of these structures 
is qualititively described and the wave propagation velocity is estimated. The problem of nonlinear 
wave stability with respect to small thermal and electromagnetic perturbations. It is shown that 
only damped perturbations correspond to space-limited solutions. 

Dissipative processes caused by Joulean heat-release in viscous motion of Abrikosov vor- 
tices give rise to different regimes of wave processes describing the final stage of evolution 
of thermomagnetic and electromagnetic perturbations in the critical state of type II super- 
conductors. A typical example of these regimes is the nonlineat stationary thermomagnetic 
wave discovered in early theoretical investigations [1-4]. It was shown in [1], that depending 
on the conditions at the supercondutor surface there may exict nonlinear stationary ther- 
momagnetic waves of two types E or H. The problem of the stability of nonlinear waves 
with respect to small thermal and electromagnetic perturbations in superconductors has not 
yet been investigated well enough. The stability of nonlinear thermomagnetic waves with 
respect to small thermal and electromagnetic perturbations in superconductors in the critical 
state is the subject of this study. It is shown that only damped perturbations correspond 
to space-limited solutions. The evolutions of the of the thermal (T) and electromagnetic 
{E, H) perturbations is described by the nonlinear equation of the thermal conductivity [1] 

z.^=A[kAT]+J^, (1) 

by Maxwel's equations 

IdH 

4:71 w 

rotH = —J (3) 
and by the equation of the critical state 

j=Ut,H)+1{E). (4) 
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where v — i^(T) and k — k{T) are the heat capacity and thermal conductivity respectively; 

— * — * 

jc is the critical current density and jV is the resistance current density. 

Let us consider a planar semi-infinite sample {x > 0) placed in an external magnetic field 
-> dH 

H — (0, 0, Hf.) growing at a constant rate = const. According to Maxwel's equation (2), 
there is a vortex electric field E = (0, E'e, 0) in the sample, directed parallel to the current 

— * — * — * 

density j: E || j; where H^. is the amplitude of the external magnetic field and E^ is the 
amplitude of the external electric field. 

Let us use the Bean-London model of the critical state for the dependence jc{T, H); — ^ = 

clH 

[5]. According to this model the dependence jc{T) can be represented as j{T) = jo[l — a{T — 
To)]: where jo is the equilibtium current density, a is the constant parameter [3] and Tq is 
the temperature of the cooling medium. 

At sufficiently high electric fields {E > Ef) the typical dependence jr{E) can be ap- 
proximated by a piecewise-linear function, jV ~ (^fE where as at small fields E < Ef the 
dependence jr{E) is, nonlinear. The nonlinearity is due to the thermo-activation creep of 
the fiow [6]. 

For the automodelling solution of the form $^ = x — vt, describing a travelling wave moving 
at a constant velocity v along the x axis, the system of equations takes the following form 

^ = (6) 



E = -H. (7) 

c 

The thermal and electrodynamic boundary conditinos for equations (5)- (7) are as follows: 

dT 

T(C ^ +oo) = To, — ^ -oo) = 0, 

^ (8) 

E{$ ^ +oo) = 0, E{$ ^ -oo) = E„ 

Excluding the variables T, H by using (5) and (7) with the boundary conditions (8) we 
obtain a differential equation for the distribution E{^): 

d^E .dE „ E^ 



Here the following dimensionless parameters are introduced z — —, [3 — L 
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Dt 

is the depth of the magnetic field penetration into the superconductor, r = — — is the 
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parameter describing the ratio of the thermal Dt 



to the magnetic D. 



coefficient; 



is the time of thermal diffusion, and 



477(7/ 

is a constant. 



diffusion 



In the approximation of weak superconductor heating (T — To) << Tq, the heat capacity v 
and the thermal conductivity k depend on temperature profile only slightly. Then at r >> 1 
we can neglect the terms with dissipative effects in equation (9) and represent its solutions 
as [7] 
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1 - th—{z - Zq) 



(10) 



Using the boundary condition E{z — > — oo) = E^. we readily determine the velocity ve'- 



E. 



2tE^ 



1/2 



of a wave with amplitude E^ and front width 

:i+t) 



5z = 16- 



-1/2 



E. 



1/2 



(11) 



(12) 



Numerical estimations give ve — i 10 



^ and 5z 
sec 



10-^ ^ 10-2 for r = 1. 



Expression (10) describes the profile of the thermomagnetic shock wave travelling into the 
supercondutor. 

To investigate the stability of the nonlinear wave with respect to small perturbations, let 
us represent the solution of the system (l)-(4) as 



T{z,t) = T{z) +5T{z,t)e^^ 



Xt 
tl 



, E{z, t) = E{z) + 5E{z, t) exp 



Xt 



(13) 



where T{z) and E{z) are the stationary solutions, and ST, 6E are small perturbations. 
Substituting (13) into (l)-(4) and assuming that ST and << T, £■ at r >> 1, which 



vv 



corresponds to a slow mstabihty A << [6] we obtain the following equations for 

determining the eigenvalues A : 



[H - A]^ = OH 



^ + ^(,),^(,)^-l + A+At/.,. 
A 



(14) 



where ^(y) = 6E{y)chy,y = —{z - zq), ^ = ^• 

The stationary wave is stable if the eigenvaluesof the operator H include negative A < 0. 
As a boundary conditions for (14) cervcs the condition that the solutions ^(y) are bounded 
at infinity. The odd term in (14) describes the effect of the thermal mode on the dynamics of 
electromagnetic perturbations in the superconductor. Using a standard procedure like that 
in [8], we can readily find the exact solution of equation (14) in the form 

l-thy\ 



*(y) = (1 - thyy-'-^ (1 + thyy+-^F a, P, 7, 



(15) 



3 



Here F is a hypergeometric function [9]; a — p + q + 3, P — p + q, ^ — 2p+l. The constant 
coefficients p and q are determined frm quadratic equations: 

It can be seen that the spectrum of cigcnfunctions of equation (14) is continuous. Analysis 
of the asymptotic behaviour of solution (15) at y ^ ±oo shows that 5E{y) is bounded only 
at A < 0. Indeed, one of partial solutions of equation (14), remaining bounded at y — > — oo 
has the form: 

5E{y) ^ eM^pVWP + q + 2,p + q-l,2p+l, exp{-2y))- 1 - 2p = uj^, (16) 

It is obvious that at ?/ +oo {exp{—2y) 0) to the asymptotics (16) corresponds a wave 
exponentially decaying at positive cOp and growing as expiiuj^y) at negative uip. It may be 
assumed that the obtained solution exhibits a regular asymptotic behaviour at y — > +oo,i.. 
the function SE{y) is bounded only at those eigenvalues for which A < 0. 

Let us represent the second linear- independent solution of equation (14) as: 

F(p + g + 2, p + g - 1, 2p + 1, s) = — -^^^^i^t^^^ X 

{p-q-l){p-q + 2) 

xF{p + q + 2,p + q-l,2p+l,l- s)+ (17) 

(2p + l)(2g) _ ^^-2,^^^ _q_i^p_q^2,-2q + l,l-s), 



{p + q + 2){p + q + l) 

where is Euler's 7 function [9]. At y — > —00 the following asymptotic representation is 
valid 

SE{y) ^ exp{ugy); 1 - 2q = coq. (18) 

It can be seen that the space-limited thermomagnetic perturbations are damped both at 
c<;p > and at Ug > 0. 

In conclusion, it should be noted that the description of the critical state used in this 
study is based on the microscopic theory of Bardeen, Cooper and Schriffer (BCS). We use 
the contunial approximation and assume that the physical characteristics of the system 
change weakly on a scale corresponding to the mean distance d between the vortices. The 
continuity condition has the form L,5z » d. A nonlinear wave can be exist in samples 
with thickness 21 satisfuing the condition Sz,L « I. 

On more restrition appears in view of the possible overheating of the superconductor to 
above the critical temperature Tg, when the system of equations (l)-(4) is not valid. Such 
a situation occurs if the condition jc(To, Eg) > a{Tc — Tq) is satisfied, Tc is the critical 
temperature. In this case, a zone of normal conductivity must appear immediately behind 
the wave front. This zone stimulates additional energy dissipation in the superconductor. 

The thermomagnetic wave E can be observed experimentally by initiating motion of a 
magnetic flux on the background of an external magnetic field growing at a constant rate 
He- Then a vortex electric field whith amplitude Eg exists in the superconductor from the 
beginning. 
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